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Abstract 



The q = 2 random cluster model is studied in the context of two mean-field 
models: The Bethe lattice and the complete graph. For these systems, the critical 
exponents that are defined in terms of finite clusters have some anomalous values 
as the critical point is approached from the high-density side which vindicates the 
results of earlier studies. In particular, the exponent j' which characterises the 
divergence of the average size of finite clusters is 1/2 and v', the exponent associ- 
ated with the length scale of finite clusters is 1/4. The full collection of exponents 
indicates an upper critical dimension of 6. The standard mean-field exponents of 
the Ising system are also present in this model (y 1 = 1/2, 7' = 1) which implies, 
in particular, the presence of two diverging length-scales. Furthermore, the finite 
cluster exponents are stable to the addition of disorder which, near the upper 
critical dimension, may have interesting implications concerning the generality of 
the disordered system/correlation length bounds. 
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A. Introduction 



The close connection between spin-systems and percolation models has lead to many impor- 
tant developments in statistical physics. For a broad class of spin models there is a mapping to 
an equivalent graphical representation having a percolation transition corresponding to the phase 
transition of the spin model (see e.g. In particular, the graphical representations of the g-state 
Potts models are the "g-state" random cluster models . For q = 2 the spin representation is the 
Ising model while the q = 1 random cluster model is ordinary bond percolation. The equivalence 
between spin models and graphical models has lead to a new class of highly efficient cluster Monte 
Carlo methods [f§] . Cluster methods simulate both the spin model and the graphical representation 
and it is often more efficient to measure thermodynamic quantities via their graphical analogs ||. 

As an example of the correspondence between spin models and graphical models, consider the 
order parameter exponent (3 that is defined for spin systems by, m ~ (T c — T)@ where T is 
the temperature, T c the critical temperature and m the order parameter. For graphical models, 
the order parameter is the fraction of sites in the percolating cluster, P^ and (3 is defined as, 
-Poo ~ (p — Pc) 13 where p is the bond occupation probability and p c the percolation threshold. For a 
given g-state Potts model and the corresponding random cluster model it turns out that m = P^ 
so that f3 = f3. Analogs of other thermal exponents may be defined for percolation models; 7 
(7') characterise the divergence of the average size of the finite cluster containing the origin below 
(above) the percolation threshold. Similarly, v (z/') characterises the divergence of average scale 
of the finite clusters below (above) threshold. The connection between the thermal and geometric 
exponents follows from the fact that the spin-spin correlation function between sites i and j is equal 
to the probability that i and j are in the same cluster in this graphical representation. Thus, in 
the high temperature phase the correlation lengths for the two models are equal and the magnetic 
susceptibility is equal to the average size of the connected clusters. Thus, in particular, 7 = 7 and 
v = v. 

On the other hand, the corresponding relation for 7' and 7' or iv' and v') is not a straightforward 
consequence of the mapping to the graphical representation. In particular, for the Potts models, 
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the graphical expression for the truncated two point function between the sites i and j consists 
of two terms: (i) The probability that the two sites are in the same finite cluster and (ii) the 
correlation of the infinite cluster density at the sites i and j. For the random cluster models (with 
q > 1) the second term is positive which leads to the inequalities v > v' and 7 > 7'. On general 
grounds, one would expect that items (i) and (ii) are comparable and for independent percolation 
(q = 1) this is the subject of a rigorous theorem |7j]. It is therefore something of a surprise that 
for percolation models of the Ising-type, the equality of 7 and 7' breaks down in mean-field. This 
was first discovered in a series of papers co-authored by one of us [P HTTll both in the context of a 
polymer/solvent model and in a mean-field calculation for the Potts/percolation droplet model. 

The principal focus of this paper is to underscore (and bolster) the conclusions inherent in the 
result 7' 7^ 7' and to clarify the relationship between the Ising magnet and its graphical representa- 
tion. In particular, the result 7' 7^ 7' in these sorts of systems will be shown in some generality: A 
full fledged Bethe lattice calculation and for the q = 2 random cluster model defined on the complete 
graph. Furthermore (at least in the context of the Bethe lattice) there is the additional mean-field 
exponent v' = 7. Taken along with (3 = (3 = ~ and 77 = 0, these exponents satisfy the standard 
scaling relations and lead to the tentative prediction d c — 6 for the upper critical dimension. Of 
particular interest is therefore the behaviour of the q = 2 random cluster model itself: a different 
set of exponents above and below threshold and the appearance of two distinct diverging length 
scales as the critical point is approached from the percolating phase. Since these models are the 
starting point of the 00 > d ^> 1 expansions, it is our belief that the behaviour uncovered here also 
holds in sufficiently high dimensions. 

It must be emphasised that the discrepancies between random cluster and Ising exponents are 
due to a difference in definitions not content. The same information is inherent in both a spin model 
and its graphical representation and it is always possible to define appropriate graphical quantities 
to represent any spin quantity. On the other hand, the fact that two natural sets of definitions 
lead to different exponents indicates the presence of real and interesting features of the Ising-type 
random cluster model in high dimensions. Some of these are discussed at the conclusion of this 
paper. 
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B. The Bethe Lattice 



Consider a half-space Bethe lattice of coordination number 3 (a binary Cayley tree). For a 
/c-level tree define the interacting bond percolation problem with the weights of configurations us 
given by W(u) = B Pe (u) q 1 ^' . Here B Pe (u) is the Bernoulli probability of to at the parameter 
value p e G [0, 1]. The quantity 1(uj) counts the number of loops. In the present context, the only 
mechanism for the formation of loops is via connections to the boundary; all points on the boundary 
are regarded as "pre-connected" f|. As is well known (for appropriate q) this is the random cluster 
representation the g-state Potts model given by the Hamiltonian —/3H = JYl<i,j>^ci,aj subject to 
the boundary condition that the boundary sites are all locked in the same state. The parameter p e 
is given by p e — [1 — e~ J ]/[l + (q — l)e~ J ] (which for the Ising case q = 2 becomes p e = tanh( J/2)). 

The partition function on the fc-level binary tree is given by Zk = XLWj^u;). There are two 
types of bond configurations: those in which the root site is connected to the boundary, and those 
in which the root site is not. Let Ik denote the total weight of configurations connected to the 
boundary and Fj, denote the total weight of bond configurations not connected to the boundary; 
Zk = Ik + Fk- Let us now take two such fc-level trees and join them by attaching their root sites to 
a new root site, thus creating a {k + l)-level tree. The recursion relations for F^+i, h+i, and Zk+i 



arc: 



Fk+i = [Fk + (1 - Pe)hf = {Z k - Pehf ; (la) 

I k+ i = 2p e I k Z k + {q- 2)p\ll ■ (lb) 

Z k+1 = Z\ + (q - l)p 2 Jl . (lc) 

We define = Ik/Zk to be the probability that a root site of a fc-level tree is connected to 

the boundary. Using the above recursion relations ([!]) and taking the limit of k — > oo, we arrive at 
the well-known expression for the percolation density: 

~ " 11111 P - " 1 + (g - l)plPl ■ [2) 
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Analysing the non-trivial solution of this equation one can see that there is a continuous phase 
transition atp e = l/2if0<g<2 with the critical exponent (3 = 1 for < q < 2 and (3 = 1/2 for 
q = 2. For q > 2 the transition becomes discontinuous. 

We now look at the average finite cluster size. Let us define a random variable c k to be the 
size of the connected component of the root site, whenever this cluster is disconnected from the 
boundary, and = if the root site is connected to the boundary. Let also X k = Y,u c kWk(u)) / Z k 
be its average. Since only the bond configurations that contribute to F k enter into this average, 
we can write X k = f k F k /Z k where f k is the average cluster size given that it is disconnected from 
the boundary. Again, merging two /c-level trees into a new (k + l)-level tree one finds the recursion 
relation for f k : 

f k+l F k+1 = (1 - Pe fZ 2 k + (1 + f k ) Pe F k (l - Pe )Z k + (1 + 2f k ) V \Fl . (3) 
Dividing both sides of (|3|) by Z k+ i and recalling the definitions of X k and we obtain 



X k+1 -1-P^ > + X k — — ^ (k) . (4) 



2^ (l-PW) + 2p e ( l-p e ) 
l + (q-l)[ Pe Pi 

Taking the limit k — > oo (the existence of which is guaranteed by a straightforward argument) we 
obtain: 

X(p e ) = lim X k = {1-PePoo? , 5) 

This expression is valid on both the ordered and disordered sides of the transition. Indeed, for 
p e < 1/2, Pqo = and Eq.@ reduces to the well known expression X = 1/(1 — 2p e ) (which gives 
7 = 1 for all < q < 2). 

The situation on the ordered side (p e > 1/2) is changed by the presence of the term proportional 
to Pqo in the denominator. For < q < 2, Poo(p e ) is asymptotically linear in (jp e — 1/2) resulting in 
the anticipated 7' = 1. However, for q = 2, 

Poo(p e ) = m(p e ) ~ (p e " ^) 1/2 » (2 Pe ~ 1) (6) 

for (p e — 1/2) small, so that the terms involving P^ are dominant. Indeed, substituting this 
expression for m directly into Eq. (||) we arrive at 



1 — /77 

X = for q = 2 and p e > 1/2. (7) 

p e m 

The result for q = 2 is 7' = 1/2 7^ 7. 

Now we look at the finite cluster distribution. The probability that the root site of a /c-level tree 
belongs to a cluster of size n that does not touch the boundary is defined as 

p(k) = F k( n ) = F k(n) ^ 
Zk F k + I k 

where F k {n) is simply the total weight corresponding to such event (clearly, F k = ^Z n F k {n)). Once 
again, merging two fc-level trees into a new one we find the following relations between these weights: 



F k+1 (l) = (l-p e ) 2 Z 2 k ; (9a) 



n-2 



F k+1 (n) = 2 Pe (l- Pe )Z k F k (n-l)+pl'£F k (i)F k (n-i-l), (n ^ 1) . (9b) 



i=i 



The next step is to divide both parts of these equations by Z k+ i and to arrive (with the help of 
(|TcD) at the expression for P^ fe+1 ^ in terms of and Pj£) . Letting k — » 00 (which is again easily 
justified) we obtain: 

Pl "l + (g-l)^' (10a) 
j n-2 

1 + (q - IjpjPi L fzi J 

It should be noted that ( |10a| ) and ( |10b| ) are independent of q when P^ = 0. Thus, below and at 
the percolation threshold (p e < 1/2) the finite cluster distribution is identical to that for the case 
of independent percolation! Since the critical exponents r (and 77) are defined at the critical point 
they must take on their mean field percolation values r = 5/2 (and r\ — 0) for any q e (0, 2]. Since 
5 = 3 for the mean field Ising model, the exponent relation r = 1/5 — 2 is violated. 

The last critical exponent of our interest here is the correlation length exponent v. In order 
to find it we adopt the standard definition of the correlation length on the Bethe lattice (see, 
e.g., Ref. 0): 



Z(Pe) = \h?r^ y Z\ x \*'dM ( n ) 



X(p 



ej x 



where T^ x (p e ) is the probability of the origin being connected to the site x while not being connected 
to the boundary. The probability Tq x depends only on the level n of the Cayley tree that the site 
x belongs to. Thus, using the metrics relation \x\ 2 = n, we can rewrite the sum in ([11]) as 



E M 2 tL = E n2 n 4 = 2 J> + l)2 n rl +1 . (12) 



|2^/ 

x 

n=0 n=0 



Performing the same procedure of merging two fc-level trees together to form a k + 1-level tree and 
taking a limit of k — > oo we generate a recursion relation for the probability r/: 

f _ f Pejl-pePpo) 

n+l ~ + ( } 

Combining this result with fll2f) we obtain the desired expression for the correlation length: 



2p e (l-p e P 0O ) 2p e X(p e ) 
\l-2p e + 2p1P 00 + (q-l)plPl h-PePoo { ] 

with X(p e ) given by Eq. ([5]). This brings us to the following relation between the critical exponents: 
v = 7/2 on both sides of the transition (in agreement with the scaling relation 7 = 2(y — rf) since 
77 = 0), which in turn leads to a surprising result: in the Ising case (q = 2) v' = 1/2 while = 1/4 ! 

So far we have dealt with the half-space Bethe lattice where the root site has only two nearest 
neighbours as opposed to three for any inner site. We claim, however, that all full-space quantities 
experience the same critical behaviour as the corresponding half-space quantities above. In fact, 
they can be explicitly calculated if we attach the root sites of two identical trees together thus 
forming a complete full-space Bethe lattice (cf. 0): 

l+p e + (g-2)p e P 00 - 
V °° ~ 00 l + (q-l) Pe Pl ' (15a) 



(15b) 



C(Pe) = £(Pe 



ej l + (g- 




3(1- 


PePoo) 


\ 2(1 +Pe 


- 2p e P 00 ) 



(15c) 



Here Poo, X, and ( refer to the isotropic, full-space quantities, which are just non-singular modifi- 
cations of the corresponding half-space quantities given by Eqs. @, (|5|) and (14). 



The disordered case. The addition of disorder has almost no effect on the previous set of results. 
This fact leads to some interesting consequences that will be discussed in the final section. Here we 
will demonstrate this stability to disorder confining attention to the exponent 7' for the case q = 2. 
The setup is as follows: The couplings are be given by (Jij) which are identical and independent 
non-negative random variables. The quantity P{£' is now a random function of these couplings that 
obeys the distributional equation 

n(^J I D ( ^ ) 

p(k+l) _ Pe,L*. oo,L T Pe,R*oc,R /, fi \ 

with P^\ and P^ R identical and independent representing the percolation probabilities for a k- 
level tree and with p Cj L and p e ,R, distributed according to tanh(J lJ /2), representing the effective 
strength of the bonds connecting the root site to the two /c-level trees situated above the root site. 
Let P^ denote the (quenched) average of P£> and similarly let p e denote the average of p e ^ or p e R . 
It is assumed that the distribution for the Jjj depends on a parameter (e.g. width, temperature) 
that can be changed continuously. The phase transition occurs at p e — 1/2 and the exponent (3 is 
the same as in the non-random case. 

Our analysis begins with the random analog of Eq. (§). After a certain amount of work, the 
relevant generalisation is seen to be 

v * u(k +l) , Pe,LX k {l -P^rPSr) +Pe,RX' k (l ~ P^lP^l) , . 

x k+1 - d i-P^ + — p{k ) p{ k) ( 17 ) 

1 + Pe,L^oo,LPe,R^oo,R 

Since the (Jjj)'s are all non-negative, all the terms in the denominator are non-negative and we 
may write 

X k+1 < d 1 - P£ + V + p e , L X k (l - Pe , R P^ R ) + p e , R X' k (l - Pe , L P£\) (18) 



Noting that all the relevant quantities in Eq. (|18|) are independent we perform the disorder average 
to obtain 

X k+1 < 1 - Pt +1) + 2p e X k (l - p e P^) (19) 

In the setup with wired boundary conditions and all the {Jij) non-negative, it is not difficult to 
show that these average quantities tend to a definite limit as k — ■> 00. Hence 

8 



Using Pqo ~ (p e — 1/2) 1 / 2 - which is not hard to prove - we obtain the (significant) first half: 
f < 1/2. 

Opposite bounds are obtained by expanding the denominator 



X k+1 > d 1 - P^ +1) + [ Pe , L X k (l - Pe , R P£> R ) + Pe, R X' k (l - p e , L P^ L )] 



X 



X [I - Pe,LPe,RP^ L P^ R \ (21) 

Neglecting positive terms and using P^lP^l — Ve^P^R — 1 we arr i ve at 

X k+1 > d 1 - Pt +1) + 2p e X k - AplP^X, (22) 

which leads to a bound similar to Eq. ( PC| ) but in the opposite direction. Hence we conclude 
f = 1/2. 



C. The Complete Graph 

Similar results can be established for the random cluster model on the complete graph. Here, 
the calculations are as straightforward as the Bethe lattice, however, the rigorous justification of 
these calculations requires some unpleasant analysis. We will again be content with the discussion 
of the exponent 7'. 

The underlying lattice consists of iV sites with bonds of uniform strength between all pairs. 
The weight for a bond configuration to is given by W(tu) = B Pe ^{uS)q l ^> oc B PN (cu)q c ^ with c 
the number of connected components and p^ defined to be 1 — e~ J l N . C.f. |13| for a more detailed 



description of the random cluster model in this context. Let G denote the size of the largest (giant) 
cluster. As is not hard to show, the probability of belonging to this cluster, G/N, converges to m(J) 
where m(J) satisfies the mean field equation m — [1 — e _Jm ]/[l + (q — l)e~ Jm ] We assume 

throughout that q = 2 and J > 2 so we are in the low-temperature phase. 

If i is a site in the graph, let Ci(N, pn) denote the analog of the quantity c k for the Bethe lattice; 
that is Ci(N,p^) is zero if i is connected to the giant cluster and otherwise is the size of the cluster 
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at i. The strategy will be to fix G and obtain estimates on (the distribution of) Cj. These estimates 
are, more or less, the desired result if the random G is replaced by mN. The large deviation result 



of in essence allows this replacement. 

Thus suppose that there are G sites in the giant cluster of an N site graph with parameter p 
(not necessarily equal to p^) and consider the cluster of the origin. For fixed e, let us assume that 
\G/N — m\ < e - otherwise for the upper bounds we will assign Cq = N and for the lower bounds, 
Co = 0. We start with the upper bounds. As before, let Fq denote the size of the cluster at given 
that it is not attached to the giant cluster. Given the condition of detachment, the origin and the 
other N — G — 1 sites act like an autonomous random cluster model subject to the condition that 
no cluster has size exceeding G - and for the upper bounds, we may neglect this condition. We will 



use the methods introduced in pi|Jl^] known as the Edwards-Sokal coupling which for complete 
graphs is particularly easy: Divide the remaining N — G sites into two groups of Aq and N 2 sites. 
One of these is identified as plus and the other as minus. There can be no bonds between spins of 
opposite type and bonds between spins of the same type occur with probability p^. Taking into 
account the relevant energetics and combinatorial factor, the result is seen to be a complete graph 
Ising problem with N — G sites and temperature parameter p = p^. Since G > (m( J) — e)N (with 
e small) it is not hard to show that the Ising system is above the critical temperature. Thus, with 
large probability, N\ is close to half of iV — G - say |2Aq — (N — G)\ < eN. Hence, when all is said 
and done, we are reduced to the problem of the cluster distribution for (subcritical) percolation on 
the complete graph with N% ~ (JV/2)(1 — m) sites and bond probability ~ J/N. Let us temporarily 
denote these parameters by n and a/n - with a < 1. As is not hard to show, the distribution for 
the cluster size at any one of these sites is bounded by a Bernoulli branching process with a mean 
of a and a maximum of n offspring. Thus Fq <d I where I satisfies the distributional equality 



/ = i + - y jtf) (23) 
n J=l 



with the IV' independent and identical in distribution to /. 

Let <7iv(e) denote the probability that \G — mN\ > eN. Let 0jv(e) denote the probability that 
|JVi - 1/2 (JV - G)\ > eN optimised over all G's such that \G - mN\ < eN. Then, solving Eq. @ 
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(after expectation) we find 



(1 - 9n)(1 - <Pn)(F ) n ,j < Ng N + (1 - 9N) 1 ^tj^ 1 m) < ^ ) 2j€ ( 24 ) 



Using the large deviation estimate of [D| for and a similar (easily derived) estimate for (p^ we 
get, letting N — > oo and then e — > 

< t _ i _ m) m 

hence 

< c » > -' s i-i^r-m) - (25b> 

The derivation of the lower bound involves a few more details. We still condition on G ~ mA^ 
and iVi ~ iV^; now we must pay lip service to the possibility of another large cluster. Let n and 
a be as before. We will imagine that there are already n — 1 sites present and that we add the 
n th at the origin. Let c satisfy c 2 n = G. Then, in order to get a cluster of size G, either (i) the 
pre-existing collection of n — 1 sites must contain a cluster of size larger than Cy/n or (ii) the new 
site must give rise to at least Cy/n bonds. For ease of future exposition, we will replace (i) by the 
weaker condition that some cluster contains at least C\fn bonds (rather than sites). Denoting by Xj 
and Xjj the indicators of these events, it is not hard to see, by comparison with the aforementioned 
branching process, that for the described values of a and n, both probabilities tend to zero at least 
as fast as exp{— b(a)y/n} for some b > 0. 

Now consider the clusters of the first n — 1 sites which we denote by K%, K 2 , . . . K s . Let ttj denote 
the probability that the new site connects to the j th cluster. Clearly p\Kj\ > ttj >p\Kj\ — p 2 \Kj\ 2 . 
Hence 

F >(1- X,)(l - X«)[l + X>l#jf -p 2 \Kj\ 3 )] (26) 

j'=i 

The first term in the sum is obviously p J2j Fj ~ essentially our upper bound. But given that Xj ^ 1, 
each \Kj\ < Cy/n and thus 



ol (y. 

F (N,p N ) > (i-X i )(l-TL u )[l + (l--j=)-'EF j (N-l,p N ). (27) 

'n n 
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The only remaining difficulty is that the F's on the right hand side of Eq. ([27D are slightly out 
of balance with regards to their arguments. However, the density p^ may be obtained from the 
density pn-i by independently removing occupied bonds with probability l/N. Thus writing Fj(N— 
IjPiv) = Fj{N — 1,ptv_i) +Fj(N — 1,Pn) — Fj(N — l,p^^i), the difference term is (distributionally) 
negative and may be bounded below by —Fj(N — l,p^_i) if even a single bond in the cluster of j 
gets removed. However since we may operate under the stipulation that there are never more than 
C\/n bonds in any of these clusters, the probability of such a loss is of the order N~ l l 2 . Putting all 
these ingredients together, we arrive at the recursive inequality 

(F } JjN >[l + ±J(l-m{J)){F )j !N }e{N,e) (28) 

with e(N,e) < 1 satisfying lim^o liniTv-^oo e(iV, e) = After a straightforward limiting argument, 
the desired result 

1 — nrn 

(Co) j = ^-rjr, T~f\\ (29) 

1 - ±J(1 - m(.J)) 

now follows from the upper and lower bounds. This gives us 7' = 1/2 for the complete graph. 



D. Conclusions/Speculations 

As indicated by our notation the exponents 7', (3 and v' have direct counterparts in spin-systems. 
A standard hyper-scaling relation, dj'/ (2/5 + 7') = 2 — rj, would therefore predict the upper critical 
dimension d c = 6 (see also Ref. |15| ) This, once again, is surprising since it differs from the usual 
value (d c = 4) associated with Ising systems. 

What has so far been demonstrated (in a tautological sense) is a breakdown in some of the 
anticipated relations between thermodynamic and geometric exponents. Let us illustrate this fur- 
ther. For magnetic systems, the critical state may be perturbed by a magnetic field h which serves 
to define the exponent 5. This exponent is related to the geometric exponent r by the following 
argument: If h <C 1, only clusters of size on the order of or exceeding 1/h will be aligned with the 
field. Hence, m(h) ~ Y,n>i/h.Pn ~ h T ~ 2 = h 1 ^ 5 , i.e. r — 2 = 1/5. This relationship has broken 



12 



down on the Bethe lattice and thus it may be presumed not to hold in sufficiently high dimension. 
It therefore must be reinterpreted as a hyperscaling relation. But there is a further point to be 
considered, namely that the above argument must also break down. 

Although this argument is far from rigorous, it appears to be irrefutable provided that one 
assumes that the distribution of clusters with h > is not significantly disturbed from the h = 
distribution and that the contribution from the infinite cluster scales as the contribution from finite 
clusters. We expect that the breakdown of r — 2 = 1/5 comes from the fact that the contribution 
from finite clusters which scale with an exponent r — 2 = 1/2 is different from the contribution of 
the infinite cluster which scales with the Ising exponent 1/5 = 1/3. 

Let us now present some speculations/conjectures concerning the behaviour of these systems in 
finite dimensions. Under the assumption that d c = 6, the exponents u', 7' etc. calculated here would 
be valid for d > 6 (perhaps with logarithmic modifications at d = 6). This would give 7^=4 = 7/3 
while Td=6 = 5/2. But what about 5 dimensions? The simplest scenario is as follows: Noting that 5 
is not a geometric exponent, let us eliminate this object in favour of rj (which is both geometric and 
thermodynamic and coincides with the Ising exponent) via the usual hyperscaling relation. This 
gives us the hyperscaling relation 

r = 3 f+ 2 -" (30) 

For d > 4 we may set 77 = 0. We note that this relation gives the right value for d — 4 and d = 6 
and breaks down for d > 6 which is consistent with the upper critical dimensionality d c = 6. We 
therefore may suppose the relation to hold for d < 6. In particular for d = 5, it gives us the 
prediction 7^=5 = 17/7. Similar reasoning leads to the predictions for 4 < d < 6 of 7' = 2/(d — 2) 
and v' = l/(d — 2), i.e. j' d=5 = 2/3 and v' d=h = 1/3. Both hyperscaling and scaling relations predict 
a! = |. On the other hand any thermodynamic definition concerning specific heat leads us to a' = 0. 
This dichotomy may be due to the fact that the singularity a! — | has a vanishing amplitude. |15 



Perhaps the most significant feature uncovered by these calculations is the appearance of two 
divergent length scales, £' and £' corresponding to the infinite and the finite clusters. Let us pause 
to reinterpret the former. Note that the (truncated) correlation function for the probability that 
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two sites belong to the infinite cluster is also the correlation function for the probability that the 
pair does not belong to the infinite cluster. Thus, £' can be related to the typical scale of cavities 
in the infinite cluster. In this light, the equality of £' and £' - as is the case for ordinary percolation 
- is eminently reasonable: Cavities in the infinite cluster of sizes up to some scale (£') that are 
populated with finite cluster that range up to a comparable scale. However our situation is quite 
different. For < T c — T <C 1 (and d > 4) we expect the cavities to be filled with relatively small 
scale finite clusters. 

Finally, let us emphasise the conclusions of the final calculations in Section B: The mean field 
exponents obtained all appear to be stable to the presence of disorder. Taken in conjunction with 
the previous discussion, - which includes the assumption/prediction that d c = 6 - this would imply 
a violation of the Harris criterion-type bound u u > 2/cf' in dimensions 5-8. Since the above 
bound holds in all systems where it is possible to define an equivalent finite-size scaling correlation 



length, |16j the implication here is that it is not possible to define such a length-scale. 

We caution the reader that the discussions in this section are highly speculative; the various 
scenarios might all be wrong. We are not definitive in any of these "predictions" - they have all been 
made without the benefit of derivations or supplementary calculations (let alone rigorous proofs). 
However in d = 5 & 6 the system under consideration is certainly within reach of currently available 
numerical methods which would shed some light on these issues. 
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